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REPORT  NO.  682 


Franz  L.  Alt 

Aberdeen  Proving  Ground,  Md. 
1  September  1948 


STEADY-STATE  SOLUTIONS  OF  THE  EQUATION  OF  BURNING 

\  ABSTRACT 

hen  a  substance  such  as  a  low -order  explosive  or  a  propellant  Is  Ignited,' 
the  burning  proceeds  with  a  speed  and  temperature  which  are  described  by  a  certain 
partial  differential  equation.  The  general  solution  of  this  equation  is  not  known. 
The  present  report  deals  with  a  special  case  of  this  equation,  namely,  the  case  of 
steady  state,  in  which  the  burning  progresses  with  uniform  speed.  We  assume  more¬ 
over  that  the  burning  substance  has  the  shape  of  a  thin  rod  of  infinite  length..  The 
limitation  to  the  steady-state  case  not  only  has  the  advantage  of  simplifying  the 
problem  mathematically,  but  it  is  interesting  because  in  many  practical  cases  the 
phenomenon  of  burning  approaches  the  steady-state  rapidly.  Thus,  the  steady-state 
solutions  presented  in  this  report  may  also  be  thought  of  as  limiting  cases  erf  more 
general  solutions.  Similarly  the  limitation  to  a  thin  infinite  rod  is  made  mostly  be¬ 
cause  it  simplifies  the  problem;  but  the  solutionsfor  many  other  shapes  do  not  differ 
much  from  the  ones  obtained  here. 


The  present  report  considers  both  the  case  of  a  perfectly  heat -insulated  rod 
and  that  of  a  rod  which  loses  heat  to  its  surroundings.  For  the  first  case,  all  pos¬ 
sible  solutions  are  obtained;  for  the  second  case,  the  report  is  limited  to  the  com¬ 
monly  observed  ranges  of  values  of  flame  speed,  heat  loss  and  room  temperatime. 
Within  these  limits,  the -various  possible  types  of  solutions  are  discussed  and  a 
small  number  of  solutions  are  presented  numerically  so  that  other  solutions  may  be 
obtained  by  interpolation,  w-—  -  >•- 

Sections  1  to  4  of  this  report  contain  the  major  results  and  an  outline  of  the 
methods  used  in  deriving  them.  Mqst  of  the  detailed  proofs  as  well  as  a  few  erf- the 
less  interesting  results  have  been  relegated  to  the  Appendix.  i ' .  ; 
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METHOD  OF  APPROACH 

It  is  desired  to  flid  the  changes  of  temperature  In  a  material  in  which  an  exothermic  chemical 
reaction,  such  as  burning,  is  talcing  place.  We  consider  the  case  in  which  the  material  has  the  shape  of  a 
thin  rod  of  infinite  length.  Each  point  on  the  rod  is  characterized  by  its  abscissa,  X,  measured  In  centi¬ 
meters  from  some  point  as  origin.  The  temperature  T  (measured  in  degrees  absolute)  in  any  point  of  the 
rod  Is  a  function  of  X  and  of  the  time  t  (measured  In  seconds  from  an  arbitrary  origin).  Unless  the  tem¬ 
perature  is  the  same  for  all  X,  beat  will  be  conducted  along  the  rod  and  will  cause  temperature  changes 
which  are  subject  to  the  well-known  equation  of  heat  flow, 

_8T  82T 

ax2 

where  k  is  a  constant  which  depends  on  the  thermic  properties  of  the  material.  At  the  same  time,  heat 

will  be  generated  In  every  point  by  the  exothermic  reaction.  We  assume  that  the  rate  of  increase  In 

-Q/RT 

temperature  brought  about  by  this  reaction  is  given  by  the  Arrhenius  expression  Ae  w  ,  where  A  and 
Q  are  constants  determined  by  the  nature  of  the  reaction,  and  R  Is  the  gas  constant.  In  addition  to  this 
rate  of  change  in  temperature  and  to  the  one  due  to  heat  conduction  along  the  rod,  we  may  have  a  certain 
loss  of  heat  to  the  space  surrounding  the  rod.  We  shall  first  consider  the  case  in  which  the  rod  Is  per¬ 
fectly  heat-insulated;  we  then  have 


^.k<^+Ae-Q/RT  (l) 

This  case  is  identical  with  that  of  a  semi -infinite  solid  hi  which  heat  flows  only  inward  from  the  surface 
and  not  parallel  to  the  surface;  and  also  with  any  other  case  of  one-dimensional  flow.  Subsequently  we 
shall  deal  with  the  case  In  which  heat  is  lost  to  the  surrounding  space  (or  gained  from  it)  by  conduction. 
In  this  case 


ST 

at 


kefT+Ae"3/RT 

ex6 


-  C  (T-Ts) 


(2) 


where  we  have  assumed  that  the  heat  conduction  between  the  rod  and  the  surrounding  medium  is  propor¬ 
tional  to  the  difference  In  temperature  (an  assumption  which  is  true  as  a  first  approximation).  C  is  the 
factor  of  proportionality,  and  T  is  the  temperature  (assumed  constant)  of  the  space  surrounding  the  burn- 
tag  rod. 


1This  assumption  implies,  somewhat  unrealistically,  that  the  burning  substance  remains  chemically  homo 
geneous  throughout  the  process  of  burning. 
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We  shall  not  attempt  to  find  all  solutions  of  the  partial  differential  equations  (1)  and  (2),  but  con¬ 
fine  ourselves  to  steady -state  solutions,  that  Is  those  solutions  for  y;hlch  T,  as  a  function  of  X  and  t,  re¬ 
mains  constant  along  all  lines 


X  -  a  t  «  const.  (3) 

r 

where  u  is  a  predetermined  constant.  In  tnis  case  the  problem  of  finding  T  for  all  X  and  t  reduces  Itself 
to  the  simpler  one  of  finding  T  for  all  X  and  one  particular  t,  say  t  -  0;  for 


T(X,  t)  -  T(X  -  ul,  0) 


We  shall  write 


(4) 


X  -  ut  >  X\  T(X  -  ut,  01  -  T*  (X*)  (4a) 

The  same  facts  may  also  be  expressed  by  saying  that  "the  points  of  constant  temperature  travel  along  the 
rod  with  uniform  velocity  u.!l 


From  (4)  we  have 


0T(X,  t)  dT*(X+)  ST(X,  t)  dT*  (X*) 

5X  dX*  '  0t  dX* 


so  that  the  partial  differential  equation  (2)  oecomes  an  ordinary  differential  equation 


dX»- 


A  e-3/Kr*  .  C  (y*  . 


(6) 


Of  the  six  parameters  in  this  equation,  three  may  be  eliminated  by  a  suitable  choice  of  the  units  of 
measurement  for  X  and  T.  If  we  set 


X*  »  X  -  ut  -  x.  VkQ/RA,  T*  (X*)  -  y(x)-  Q/R 

u  -  v  KkAR/Q,  C  -  c-AR/Q,  Tg  -  yg  •  Q/R 


(5a) 


equation  (5)  becomes 


i  «■  v  g  +  e-Vy  -  c  (y  -  y  )  -  0  .  (6) 

dx 

The  symbols  y  and  y  represent  absolute  temperatures,  expressed  In  units  of  Q/R.  Similarly,  v  is  the  ve¬ 
locity  of  propagation  of  the  burning  process,  and  x  may  be  considered  as  either  time  or  distance,  all  meas¬ 
ured  in  suitable  units. 
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Equation  (6)  Is  the  "equation  of  Bteady-atate  burning."  In  order  to  solve  It  we  note  that  It  does  not 

contain  x  explicitly.  Therefore,  Us  order  can  be  reduced  by  one.  We  consider  y,  rather  than  x,  as  the  in- 

2  2 

dependent  variable  of  the  problem  and  set  dy/dx  -  -z(y).  Then  dz/dy  *  (dz/dx)  •  (dx/dy)  -  (d  y/dx  )  •  (l/o), 
and  (8)  becomes 

»  v  .  I  (e_1/y  -  c  (y  -  ys) )  (7) 

This  is  a  non-linear  ordinary  differential  equation  oi  the  first  order.  We  shall  refer  to  It  as  the  "trans¬ 
formed  equation."  Our  principal  task  will  be  to  solve  it;  this  Is  accomplished 'by  numerical  Integration. 
Once  this  Is  done,  the  solution  of  (6)  can  be  reduced  to  a  simple  quadrature,  as  follows:  instead  of  finding 
y  as  a  function  of  x,  we  find  x  as  a  function  of  y,  from  the  equation 


‘  j'S) 


where  z(y)  is  the  solution  of  (7).  Having  solved  (6),  we  obtain  the  solution  of  (5)  by  changing  the  scale  factors 
as  in  (oa),  and  the  solution  of  (2)  by  means  of  (4). 

Equations  (5)  and  (6)  are  of  the  second  order,  and  therefore  two  boundary  conditions  are  required  in 
order  to  determine  a  particular  integral  of  either  of  these  equations.  One  of  these,  which  may  be  given  in  the 
form  of  prescribing  the  values  of  the  temperature  or  Its  first  derivative  for  a  particular  value  of  x,  only 
serves  to  fix  the  location  of  the  solution  with  respect  to  the  x-axis.  This  is  because,  If  one  solution  is  given. 
Infinitely  many  others  may  be  obtained  by  shifting  the  former  by  an  arbitrary  amount  parallel  to  the  x-axls. 
The  olhe^  boundary  condition,  however,  has  an  essential  Influence  upon  the  shape  of  the  solution  curve.  It  Is 
customary.  In  dealing  with  the  equation  of  burning,  to  prescribe  that  the  process  of  burning  should  start  from 
a  non-burning  equilibrium  state:  In  the  case  without  heat  less,  the  burning  substance  Is  supposed  to  have 
originally  been  at  zero  temperature  (since  at  any  other  temperature  tt  is  not  in  equilibrium) ,  and  in  the  case 
with  heat  loss,  it  is  supposed  to  have  originally  been  at  that  temperature  y  which  is  just  so  much  higher 


than  the  room  temperature  yg  that  the  heat  generated  within  the  substance,  e 


1  o,  Is  balanced  by  the  heat 


lost  to  the  surrounding  space,  c(yQ  -  y  ).  In  other  words,  it  is  customary  to  prescribe,  as  the  second 
boundary  condition,  that  y  -  0  (in  ‘.he  ease  e  ®  0)  or  y  -  (in  the  ease  c/0)  "a  long  time  ago."  It  Is  not  pos¬ 
sible  to  postulate  these  conditions  for  any  finite  time  (nor,  by  the  same  token,  for  any  finite  value  of  x), 
since  this  would  contradict  the  assumption  of  a  steady  state.  Therefore,  the  second  boundary  condition  Is 
formulated  by  stating  that  the  temperature  y  should  approach  the  limit  stated  (0  or  y  )  as  x  tends  to  Infinity, 
or,  which  Is  the  same,  as  t  tends  to  -  bo  .  Evidently,  z  -  -dy/dx  must  simultaneously  tend  to  zero.  Thus, 

In  terms  of  the  transformed  equation,  this  boundary  condition  states  that  z  »  0  for  y  «  0  (in  the  case  c  «  0) 
or  for  y  -  yQ  (in  the  case  c  /  0). 


SOLUTIONS  or  THE  TRANSFORMED  EQUATION  WITHOUT  HEAT  LOSS 

We  consider  first  the  case  c  -  0,  i.e..  the  case  of  a  perfectly  insulated  rod.  We  confine  ourselves 
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to  non-negative  values  of  y.  This  limitation  suggests  itself  because  of  the  Interpretation  of  y  as  the  absolute 
temperature  of  a  substance  (except  for  a  scale  factor).  Figure  1  shows  the  shape  of  several  solution  curves 
of  (7).  Figures  1-4  are  not  drawn  to  scale.  They  merely  show  the  general  appearance  of  the  solutions. 

In  Figures  1,  2  and  4  the  temperature  (y)  is  plotted  horizontally,  and  the  rate  of  change  of  temperature  (z)  is 
plotted  vertically. 

There  are  infinitely  many  solution  curves  similar  to  curve  I-A  of  the  figure.  Tney  start  at  the  origin 
with  horizontal  tangent,  increase  to  a  maximum,  then  decrease,  intersect  the  y  axis  vertically,  turn  back  In 
the  direction  of  decreasing  y  and  reach  the  z  axis  for  a  certain  negative  value  of  z;  at  the  intersection  with 
the  z  axis  the  slope  is  v.  There  is  one  and  only  one  curve  of  this  type  through  every  point  in  the  area  bound¬ 
ed  by  the  a  axis  and  curve  I-B. 

Cuive  I-B  is  similar  in  shape  to  the  curves  of  Type  I-A  except  that  Its  slope  at  the  origin  is  v 
rather  than  0.  There  is  only  one  curve  of  this  type.  In  the  application  of  this  differential  equation  to  the 
problem  of  burning  this  Is  the  most  interesting  solution. 

The  solutions  of  Type  I-C  are  similar  in  shape  to  curve  I-B  but  they  start  at  a  point  on  the  z  axis 
with  some  positive  value  of  z,  rather  than  at  the  origin.  The  slope  at  this  point  is  again  v.  There  are  in¬ 
finitely  many  solutions  of  this  type  One  and  only  one  passes  through  each  point  in  the  area  between  curves 

I- B  and  n-C. 

Solution  It-C  starts  at  the  z  axis  with  positive  z  and  slope  v,  increases  monotonically  and  approaches 
the  value  z  -  1/v  asymptotically.  There  Is  only  one  solution  of  this  type. 

The  solutions  of  Type  in -C  start  at  the  z  axis  with  a  value  of  z  larger  than  that  of  solution  II-C. 

Their  slope  at  the  starting  point  is  v.  With  Increasing  y  they  Increase  to  infinity.  Their  slope  decreases 
from  y  to  a  positive  minimum  value,  then  Increases  again  and  approaches  v  as  y  goes  to  infinity.  One  and 
only  one  solution  of  this  type  passes  through  each  point  in  the  area  above  curve  II-B.  The  curves  of  Type  I-C, 

II- C  and  III-C  satisfy  the  differential  equation  but  not  the  boundary  condition  z(0)  »  0.  On  the  other  hand,  all 
curves  of  Type  I-A  as  well  at  the  single  curve  of  Type  I-B  are  admissible  solutions  of  our  problem. 

The  five  types  of  solution  curves  which  we  have  described  and  which  are  sketched  in  Figure  1  exist 
if  v  Is  considerably  smaller  than  one.  If  the  value  of  v  Is  increased,  the  starting  point  of  solution  n-C'  moves 
closer  to  the  origin.  For  a  certain  value  of  v  It  coincides  with  the  origin.  For  this  v  there  are  no  so.utlons 
of  Types  I-B  and  I-C.  We  shall  denote  this  limiting  velocity  by  V.  Its  value  has  been  determined  by  re¬ 
peated  numerical  integration: 

V  *  .90280. 

If  v  increases  beyond  this  value,  then  the  solution  curves  have  the  shapes  shewn  in  Figure  2.  The 
solutions  of  Type  I-A  are  similar  to  those  for  small  v.  They  fill  the  area  below  curve  II-A.  The  latter  is 
a  solution  curve  which  starts  at  the  origin  with  slope  0  and  increases  asymptotically  toward  z  »  1/v. 

Curves  U3-A,  III-B  and  III-C  of  Figure  2  are  similar  to  the  curves  ni-C  of  Figure  1,  except  that  TTI-A  and 

III- B  start  at  the  origin  and  m-A  has  a  horizontal  slope  at  the  origin.  There  is  only  one  solution  of  Type 
III-B  but  infinitely  many  of  Types  in-A  and  III-C. 
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SOLUTIONS  OF  THE  EQUATION  OF  STEADY-STATE  BURNING  WITHOUT  HEAT  LOBS 

In  the  preceding  paragraph  we  have  discussed  solutions  of  the  transformed  Equation  (7).  We  now 
turn  to  the  corresponding  soluTcns  of  Equation  (G)  still  limiting  ourselves  to  the  case  of  perfect  Insulation 
(c  -  0).  These  solutions  can  be  found  by  the  simple  integration  indicated  in  Equation  (8). 

Solutions  of  Equation  (6)  corresponding  to  the  five  types  of  Figure  1  are  shown  In  Figure  3.  Note 
that  on  this  diagram  the  temperature  y  is  plotted  vertically  whereas  in  Figure  1  it  was  plotted  horizontally. 

The  horizontal  axis  in  Figure  3  may  be  Interpreted  in  two  ways:  as  the  X-axis  along  the  burning 
rod,  where  burning  is  assumed  to  progress  from  left  to  right;  or  sis  the  time  axis,  where  time  is  read  from 
right  to  left.  In  the  termer  Interpretation,  each  of  the  curves  shown  in  Figure  3  is  a  possible  temperature 
distribution  along  the  rod,  "possible"  meaning  a  distribution  which  will  result  in  a  steady  flow  of  heat;  that 
is  to  say,  if  one  of  these  distributions  exists  initially,  then  the  generation  of  heat  within  the  material  of  the 
rod  and  the  flow  of  neat  along  the  rod  will  change  the  temperature  in  each  point  in  such  a  way  that  after  any 
Interval  of  time  the  temperature  distribution  will  appear  to  have  been  shifted  to  the  right  without  having 
changed  its  shape.  In  the  second  Interpretation,  we  may  find  out  how  the  temperature  in  any  desired  point 
of  the  rod  changes  with  time,  by  traveling  along  one  of  the  curves  from  right  to  left.  We  see  that  "a  long 
time  ago"  the  temperature  was  zero  (absolute),  that  it  rises  gradually,  and  that  on  some  curves  it  reaches 
a  maximum  and  declines  again. 

The  solutions  of  Types  I-C,  II-C  and  III-C  reach  an  absolute  temperature  of  0  at  some  finite  point 
X  (or,  in  the  other  interpretation,  the  temperature  starts  at  0  at  some  finite  time).  These  solutions  are, 
therefore,  not  representations  of  any  steady-state  phenomenon.  Solutions  of  Type  I-A  or  I-B,  on  the  other 
hand,  can  be  interpreted  as  temperature  curves  along  a  burning  rod  in  the  steady -state.  Oh  each  of  these 
curves  y  has  a  maximum  which  might  be  thought  of  as  the  temperature  of  the  burning  front.1  This  front 
should  be  visualized  as  proceeding  along  the  burning  rod  from  left  to  right  (the  direction  of  burning  is  de¬ 
termined  by  the  sign  of  u  In  Equation  (3) ).  The  portion  of  the  temperature  curve  to  the  left  of  the  burning 
front  (and  correspondingly  that  portion  of  the  solutions  of  (7)  for  which  z  Is  negative)  covers  that  part  of 
the  rod  which  Is  already  burned.  Because  of  the  chemical  changes  caused  by  burning,  Equation  (1)  no 
longer  represents  the  temperature  changes  m  this  region. 

in  Figure  3  all  solution  curves  have  been  drawn  through  one  point.  That  is  to  say,  we  assume 
that  the  temperature  In  that  particular  point  of  the  rod  is  given.  Other  solutions  are  obtained  by  shifting 
those  shown  In  the  diagram  to  the  left  or  right  by  any  desired  distance. 

Certain  points  in  Figure  3  are  lettered  the  same  as  the  corresponding  points  In  Figure  1,  to 
facilitate  comparison. 

The  maximum  ordinate  of  curve  i-B  is  of  particular  interest.  It  represents  the  highest  burning 


Actually  this  temperature  will  usually  differ  from  the  observed  burning  temperature  because  the  conditions 
prevailing  in  a  turning  substance  are  not  accurately  described  by  Equation  (1),  especially  not  at  high  tem¬ 
peratures  at  which  most  burning  substances  undergo  chemical  changes.  It  Is  probable  that  a  point  somewhat 
to  the  right  of  the  maximum  in  Figure  3  corresponds  more  closely  to  the  actual  burning  front. 
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temperature  which  could  possibly  be  reached  In  steady-state  burning.  Its  value  depends  on  v.  The  following 
values  were  determined  by  numerical  integration: 

v  Maximum  y 


10-1 

0.2 

►—* 

O 

1 

CO 

0.0692 

-4 

2  x  10* 

0.0565 

2  x  ID*5 

0.044  8 

2  x  10‘8 

0.0277 

These  values  are  plotted  in  Figure  5.  The  accuracy  of  the  last  decimal  place  is  doubtful.  The  figures  can 
be  converted  to  customary  units,  (degrees,  centimeters,  etc.)  by  means  of  Equations  (5a).  As  v  gro.vs  be¬ 
yond  0.1  toward  the  value  V  »  .90280,  the  maximum  of  y  grows  rapidly  to  infinity. 

The  only  difficulty  in  the  integration  process  (8)  arises  when  z  becomes  0.  This  occurs  at  most 
twice  on  each  solution  curve.  It  can  be  shown  lhat  at  the  points  where  the  solution  of  (7)  intersects  the  y 
axis  with  vertical  slope,  the  integral  of  1/zhas  a  finite  value.  On  the  other  hand,  for  those  solutions  of  (7) 
which  pass  through  the  origin,  the  i.Tegral  of  1/z  gees  to  infinity  as  y  tends  to  0. 

SOLUTIONS  ALLOWING  FOR  HEAT  LOSS 

We  now  turn  to  the  case  in  which  c  has  some  positive  value,  that  is  the  case  In  which  heat  is  lost 

to  the  medium  surrounding  the  burning  substance.  The  shapes  of  the  solution  curves  In  this  case  vary 

greatly  with  the  varying  values  of  v,  c  and  yg.  In  this  discussion  we  shall  confine  ourselves  to  those  values 

which  are  found  in  actual  experiments  with  burning  substances.  The  range  for  v  is  that  shown  in  the  tabie 

*15  -22 

in  the  preceding  section.  Realistic  values  of  c  are  of  a  magnitude  between  10  and  10  .  For  y  ,  which 

represents  room  temperature,  the  values  0.012  and  0.017  are  used.  The  difference  between  these  values  is 
caused  not  so  much  by  variations  in  real  room  temperature,  as  by  the  differences  in  the  constant  Q,  which 
determines  the  scale  in  which  y  is  measured. 

Furthermore,  this  discussion  will  be  confined  to  that  range  of  y  which  corresponds  to  temperatures 
commonly  measured.  An  added  reason  for  this  limitation  i3  that  for  extreme  '’clues  of  y  the  heat  loss  Is 
probably  not  adequately  represented  by  the  term  c(y  -  y  ).  Some  Indication  of  the  behavior  of  solutions  in 
this  extreme  range  is  given  in  the  Appendix. 

The  shape  of  the  solution  curves  of  (7)  is  shown  in  Figure  4,  In  this  figure  yQ  represents  the  tem¬ 
perature  of  the  rod  before  burning,  y1  the  Ignition  temperature.  The  vslus  yQ  is  slightly  larger  than  room 
temperature.  This  is  because  even  at  room  temperature  some  heat  is  generated  in  the  combustible  sub¬ 
stance  by  the  chemical  reactions  taking  place.  The  difference  in  temperature  between  the  rod  and  the  sur¬ 
rounding  space  must  be  just  large  enough  so  that  the  small  amount  of  heat  generated  equals  the  heat  lost 
to  the  surrounding  space.  As  the  temperature  rises  above  y  ,  the  heat  loss  increases  at  first  more  rapidly 
than  the  amount  of  heat  generated.  For  still  higher  temperatures,  however,  the  latter  increases  more 
rapidly;  when  the  temperature  reaches  the  value  y^,  heat  loss  and  heat  generation  are  again  in  balance. 
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As  before,  we  shall  classify  the  solutions  according  to  their  behavior  with  increasing  y.  Some 
reach  a  maximum  and  then  decrease  to  z  *  0,  (Type  I),  some  grow  to  infinity  (Type  III),  and  one  approaches 
asymptotically  the  value  z  -  1/v  (Type  II).  The  further  subdivision  according  to  shape  near  the  singular 
points  differs  somewhat  from  the  case  without  heat  loss. 

In  Figure  4,  to  the  right  of  y^  the  shape  of  the  curves  is  similar  to  those  in  Figure  1.  However, 
the  slope  of  Type  I-A  curves  at  the  point  y  Is  not  0  but  a  small  positive  constant.  The  slope  of  curve  I-B 
at  the  point  y1  is  slightly  smaller  than  v.  All  other  curves  Intersect  the  line  y  =  y^  with  slope  v. 

To  the  left  of  there  are  infinitely  many  solution  curves  through  the  point  y  »  y^,  z  »  0,  All  of 
them  have  the  same  small  positive  slope,  except  for  one  which  has  a  slope  slightly  smaller  than  v.  The 
latter  is  the  continuation  of  curve  I-B,  the  former  are  the  continuations  of  the  curves  of  Type  I-A.  Of 
these,  one  passes  through  the  point  y  -  y  ,  z  -  0,  while  the  others  are  above  or  below  It. 

In  order  to  have  our  boundary  conditions  satisfied,  we  have  to  consider  solutions  in  the  (y,  z) 
plane  through  the  point  y  -  y  ,  z  =  0.  There  are  only  two  solution  curves  through  this  point,  one  increasing 
■and  one  decreasing.  This  fact  is  of  considerable  importance  in  the  physical  interpretation  of  the  solutions. 
The  decreasing  curve  represents  not  a  burning  substance  but  rather  a  substance  in  which  the  temperature 
is  initially  at  the  ignition  point  and  drops  slowly  toward  room  temperature.  (In  Figure  4  this  is  the  curve 
of  Type  I-A  mentioned  above;  with  a  different  selection  of  parameters  it  might  be  a  curve  of  different 
type.)  T.ie  other  of  the  two  curves  passing  through  the  point  yn  is  marked  "I-D"  in  Figure  4.  It  represents 
steady-state  owning.  It  is  Interesting  to  note  that  the  latter  solution  is  unique.  On  the  other  hana,  in  the 
case  without  heat,  loss  (c  -  0)  we  saw  that  there  are  Infinitely  many  solutions  which  start  at  0  temperature 
and  represent  steady-state  burning. 

T  he  Infinitely  many  solution  curves  passing  through  the  point  y^  represent  substances  in  which 
the  temperatwe  was  initially  at  the  ignition  point.  Since  this  is  a  point  of  unstable  equilibrium,  it  is  not 
surprising  to  find  that  there  are  infinitely  many  ways  In  which  the  temperature,  starting  from  this  value, 
can  either  increase  or  decrease. 

Flgwe  4  shows  the  field  of  the  transformed  equation  (7),  similar  to  Figures  1  and  2.  It  Is  easy  to 
visualize  the  shape  of  the  corresponding  solutions  of  the  equation  of  burning,  in  analogy  to  Figure  3.  In 
particular,  the  solution  corresponding  to  I-D  in  Figure  4 -the  only  one  which  satisfies  our  boundary  con- 
ditions— looks  like  I-B  in  Figure  3,  except  that  its  asymptote  on  the  right  is  not  zero  but  yQ. 

It  Is  of  primary  interest  to  determine  numerically  the  maximum  temperatures  for  the  unique  so¬ 
lution  of  steady-state  bwning.  The  following  values,  were  obtained  by  numerical  integration  of  Equation  (71 


V 

ys 

Maximum  y 

10-3 

.012 

,0672 

10*3 

.017 

.0663 

2  x  10'5 

.017 

.0428 

2  x  10'8 

.012 

.0288 

These  figures,  as  well  as  those  given  in  "Solutions  of  the  Equation  of  Steady -State  Burning  without  Heat  Loss,' 
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previously,  are  plotted  In  Figure  5.  The  accuracy  of  the  last  decimal  place  is  In  doubt.  The  difference  be¬ 
tween  the  values  of  y  given  in  the  two  tables  is  approximately  3%  in  the  case  where  yg  is  .012,  and  4%  where 
yg  is  ,017.  These  percentages  can  certainly  not  remain  constant  for  all  values  of  v,  but  in  the  range  covered 
by  this  report  they  represent  a  sufficiently  good  approximation. 

Note  that  the  value  of  c  is  not  listed  in  the  above  table.  This  is  because  the  results  are  practically 
Independent  of  c  as  long  as  c  is  small  but  positive.  The  presence  of  c  does,  however,  have  the  measurable 
effect  of  "stopping'1  the  solution  curves  at  y  -  y  ,  rather  than  at  y  -  0. 

The  temperature  curves  which  have  been  calculated  numerically  are  shown  In  the  attached  graphs 
9-13.  Those  for  c  =  0  are  solutions  of  Type I-B,  those  for  positive  e  are  solutions  of  Type  I-D.  Figures  9-12 
are  plots  of  the  temperature  gradient  against  temperature,  similar  to  Figures  1  and  4.  Figure  13  shows 
the  temperature  distribution  along  the  burning  substance,  similar  to  Figure  3.  One  difference  between 
Figures  3  and  13  Is  that  the  former  shows  all  integral  curves  going  through  one  point,  whereas  in  the  latter 
each  curve  is  assumed  to  reach  its  maximum  for  x  »  0.  This  is,  of  course,  a  matter  of  choice  of  initial 
conditions.  Moreover,  Figure  13  has  a  logarithmic  scale  for  x,  so  that  x  *  0  is  not  shown. 
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APPENDIX 

Computational  Procedure.  In  the  application  oi  Equations  (1)  and  (2)  to  steady-state  burning,  the 
most  Interesting  solutions  are  (a)  in  the  case  c  »  0,  the  solution  of  Type  I-B,  and  (b)  In  the  case  of  positive 
c,  the  solution  of  Type  I-D,  In  the  form  (7)  the  equation  can  be  integrated  numerically  by  one  of  the  usual 
methods.  Only  the  points  for  which  z  =  0  need  special  attention,  since  in  these  points  the  right-hand  side  of 
(7)  is  not  determined.  These  points  shall  now  be  considered. 

We  shall  use  the  abbreviations 


p(y) 


-i/y 


-  c  (y  -  ys) 


(9) 


q  (y,  z)  =  v 


p(y) 

z 


(10) 


so  that  Equation  {'/)  can  be  written 

'g-q(y,z>  (id 

The  points  in  which  the  integrals  of  this  equation  intersect  the  y-axis  are  of  two  types:  those  where  p(y)  /  0, 
and  those  where  p(y)  =  0. 

In  the  former  case,  if  we  approach  a  point  z  =  0  along  an  integral  curve  of  (11),  dz/dy  increases  be¬ 
yond  all  bounds.  The  process  of  numerical  integration  can  nevertheless  be  continued  in  the  neighborhood  of 
these  points  by  changing  variables,  considering  z  as  the  independent  and  y  as  the  dependent  variable,  and  in¬ 
tegrating  the  equation  in  the  form 


dy  _ _ 1 _ 

dz  1  ,  -l/y  ,  \ , 

v  -  -  (e  3  -  c  (y  -  yg) ) 

It  remains  to  Investigate  the  points  where  z  »  0  and  p(y)  -  0.  it  can  be  shown  by  elementary  con¬ 
siderations  that  the  function  p(y)  vanishes  in  one,  two  or  three  points  (for  real  positive  y)  depending  on  the 

values  cf  the  constants  c  and  y  .  If  c  is  larger  than  or  equal  to  4  e  ,  then  p  vanishes  In  one  place  only, 

s  -2 

regardless  oi  the  value  of  y  .  If  c  is  less  than  4  e  ,  then  p  vanishes  in  three  places  provided  y  lies  be- 

s  -1  s 
tween  certain  boundaries.  For  c  less  than  e  ,  the  lower  of  these  boundaries  is  zero,  i.e.,  p  vanishes  in 

tiiree  place?;  provided,  only  that  y  is  small  enough.  In  practice,  c  is  always  very  much  smaller  than  either 

s  -15  -22 

of  the  two  limits  mentioned,  being  of  a  magnitude  between  10  and  10  .  Moreover,  in  all  practical 

cases  ys  is  so  small  that  p(y)  vanishes  for  ttiree  different  values  of  y.  We  shall  denote  these  three  places 
in  order  by  y^,  y^,  y^.  The  smallest,  yQ,  Is  always  slightly  larger  than  yg;  it  is  the  temperature  of  the 
combustible  material  before  burning  has  begun,  y^  is  the  ignition  temperature,  y^  is,  in  the  usual  appli¬ 
cations,  so  large  that  it  is  never  reached. 

L:  c  =  0,  then  p(y)  vanishes  only  for  y  ■*  0. 

We  shall  see  below  that  if  an  integral  curve  of  (11)  passes  through  the  point  (y^ ,  0),  then  in  that  point 


BBJF' 
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Its  slope  has  one  of  the  two  values 


1 


If  an  Integral  curve  of  (11)  passes  through  the  point  (y  ,  0),  then  its  slope  in  that  point  has  one  of  the  two 
values 


In  practical  cases,  these  two  values  are  almost  equal  to  v  +  ^  and  to  -  p  respectively.  If,  in  the  case  c  *  0, 
an  integral  curve  passes  through  the  origin,  then  its  slope  at  the  origin  has  one  of  the  two  values  v  or  0. 

For  the  numerical  computation  It  is  important  to  remember  that  the  integral  curve  going  through 
the  origin  with  slope  v,  or  through  (yQ,  01  with  slope  v  +  c/v,  remains  almost  straight  for  a  considerable 
distance.  Thus,  the  numerical  Integration  along  these  curves  can  be  begun  with  fairly  large  steps.  At  some 
distance  from  the  origin  or  the  point  (yfl,  0),  the  slope  begins  to  differ  noticeably  from  its  initial  value,  and 
then  the  length  of  the  integration  step  must  be  diminished. 

The  solution  of  Equation  (8)  is  obtained,  as  explained  in  "Method  of  Approach"  (page  8),  by  in¬ 
tegrating  the  differential  -  .  However,  the  integrand  becomes  infinite  when  2  is  zero,  if  simultaneously 

p(y)  is  different  from  zero,  then  the  Integration  can  nevertheless  be  carried  out;  for  in  the  neighborhood  of 
such  a  singularity 


dy  11.  dz  dz 

z  '  z  '  dz  '  ,  Elzi)  *  z  v  -  p(y'> 

dy  '  z  ’ 

and  in  the  last  form  the  denominator  differs  from  zero  even  at  the  singular  point.  On  the  other  hand,  if  p(y) 
also  vanishes,  then  It  can  easily  be  seen  that  the  Integral  becomes  infinite.  In  other  words,  corresponding 
to  those  solutions  of  (7)  which  go  through  the  origin  or  through  the  point  (yo,  C)  or  (y^  0),  we  have  solutions 
of  (6)  which  approach  asymptotically,  as  x  goes  to  infinity,  the  value  y  -  0  or  y  -  yo  or  y  -  y1> 

Integral  curves  in  the  singular  points  (z  =  p(y)  =  0).  (a)  In  the  case  c  /  0,  the  behavior  of  the  inte¬ 
gral  curves  of  (7)  or  (11)  can  be  determined  by  known  methods  (see,  for  example,  E.  Goursat,  Course  in 
Mathematical  Analysis,  Vol.  II,  Part  II,  pp.  179-180).  If  (y*,  0)  represents  one  of  the  three  singular  points 
(y  >  °)>  (y,,  0),  (y„,  0)  of  Equation  (11),  the  Integral  curves  of  (11)  in  the  neighborhood  of  this  point  are 
known  to  be  similar  to  those  of  the  equation 


dz  _  v  z  -  (y  -  y»?  ■  p‘  (y*) 
3y  z 

Integral  curves  of  the  latter  equation  have  at  (y*,  0)  one  of  the  two  slopes 

s'  *  |v  +  y'iv2  -  p*  (y*) 


14 


BALLISTIC  RESEARCH  LABORATORIES 


and 


s"  -  |v  -  |/jva  -  p*  (y‘) 

If  s',  s"  are  real  and  have  opposite  sign  (i.e.,  if  p'  Is  negative,  which  is  the  case  for  and  but  not  for  y^), 
then  there  are  only  two  integral  carves  passing  through  the  singular  point;  the  latter  is  called  a  "Saddlebadc". 
(See  the  point  yo  in  Figure  4).  At  (y^,  0),  p'  is  positive,  so  that  s',  s"  are  either  real  and  positive  or  imagi¬ 
nary.  If  they  are  real,  the  point  is  a  node;  through  it  pass  one  integral  curve  with  slope  sf  and  infinitely  many 
with  slope  s".  If  s',  s"  are  complex  (or  real  but  equal),  then  the  point  (y1 ,  0)  Is  a  focus,  the  Integral  curves 
have  the  form  of  spirals  which  approach  the  focus.  The  corresponding  solutions  of  equation  (6)  give  a  tem¬ 
perature  which  oscillates  above  and  below  the  Ignition  point  y^.  The  condition  for  the  existence  of  such  os¬ 
cillating  solutions  is 

p*  (yt>  > ■ 

(b)  In  the  case  e  =  0  the  above  method  is  not  applicable  because  p*  vanishes  for  y  =  0  (which  is  the 
only  singular  point  In  this  case).  We  can,  however,  proceed  as  follows: 

Suppose  that  a  solution  of  (7)  or  (11)  goes  through  the  origin,  and  has  in  that  point  a  derivative 
z'  (0).  Then 

z*  (0)  =  lim  5.M 
y=0  y 

According  to  a  theorem  about  real  functions,  there  Is  a  sequence  of  points  y,  converging  to  0,  for  which 


z*  (0)  =  lim  z*  (y)  . 
y-0 


Thus 


z'  (0)  =  lim  (v  -JLiZ)) 
y=0  Z 


v  -  lim  (P  <■  v  -  lim  -E-lZ]  .  lim  ^ 

y-0  y  Z  y =0  ^  y*0  Z 


(12) 


The  last  step,  however,  is  permissible  only  if  we  know  that  at  least  one  of  the  two  limits  on  the  right  con¬ 
verges  to  a  value  different  from  zero.  Now  if  z'  (0)  /  0,  then 


lim  £ 
n  Z 
y=0 


exists  and  is  equal  to  1/z1  (0),  and  therefore  different  from  0,  and  it  follows  that  the  last  equation  of  (12) 
holds.  Its  only  solution  is  z‘  (0)  »  v.  Thus,  z'  (0)  =■  0  and  z‘  (0)  -  v  are  the  only  possible  slopes  of  integral 
curves  of  (7)  at  the  origin  in  the  case  c  *  0. 
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We  shall  prove  In  the  next  section  that  the  origin  is  a  node;  there  is  only  one  integral  curve  of  (?) 
going  through  the  origin  with  slope  v,  while  there  are  infinitely  many  with  slope  0.  This  proof  is  based  on  a 
study  of  the  curvature  of  the  integral  curves. 

The  Curvature  of  Integral  Curves.  In  the  further  discussion  we  .-hall  denote,  as  before, 


P'  M  ‘  ^ 


(13) 


and  furthermore  by  q'  the  derivative  of  q  with  respect  toy  in  the  direction  of  the  integral  curve  of  (7); 


a*  (V  2) „  gq  (y. 2'  +  gqJy.  ^ 

q  ly’  •  dy  +  6z  dy 


q  (y.  t) 


(141 


An  integral  curve  going  through  the  point  (y,  z)  Is  convex  or  concave  in  that  point  depending  on  the  sign  of 
q‘  (y,  z).  In  the  (y,  z)  plane  the  region  where  q*  is  positive  is  separated  from  the  region  where  it  is  negative 
by  a  curve  along  which  q*  -  0  (and  possibly  by  portions  of  the  y-axls,  on  which  q1  is  not  defined). 

For  this  curve  we  have 


0»q' 


P 

S 

z 


or 


vp-  z  +  p 


(15) 


This  equation,  which  determines  z  as  a  function  of  y,  Is  the  equation  of  the  curve  on  which  q’  -  0.  The  dis¬ 
criminant  of  this  quadratic  equation. 


D  «  p^  (v2  -  4p5) 

2 

vanishes  if  either  p  -  0  or  v  =  4p\  It  is  Important  to  find  the  places  where  D  »  0  because  Equation  (15)  gives 

«) 

two,  one  or  no  value  of  z  for  a  given  y  depending  on  the  sign  of  D  for  that  y.  Now  p  vanishes  in  the  places 
70,  y-,  and  y2  in  the  case  c  /  0,  and  for  y  -  0  only  In  the  case  c  >  0;  everywhere  else  p2  is  positive.  On  the 
other  hand,  the  sign  of  vS  -  4p'  depends  on  the  value  of  v  and  c.  The  function  p*  (y)  has  a  maximum  at 
y  -  1/2;  Its  value  Is  p'  (j)  »  4  e  -  c.  From  this  value,  p'  decreases  to  both  sides,  vanishes  in  two  places 
which  we  shall  call  y*  and  yg,  and  approaches  the  value  -c  as  y  approaches  0  or  infinity.  Thus,  if  v  is 
larger  than 


-2 
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p 

then  v  -  4p‘  is  always  positive,  and  D  is  positive  everywhere  except  where  p  =  0.  In  this  case  the  curve 
q'  -  0  consists  of  two  branches  which  are  separate  for  all  positive  values  of  y,  except  y  ,  y^,  y^.  In  the  two 
places  where  p'  ■*  0,  Equation  (15)  becomes  linear  and  has  only  one  solution,  namely,  z  »  p/v;  one  of  the  two 
branches  of  the  curve  q*  »  0  goes  to  z  =  co .  The  shape  of  the  curve  is  as  shown  in  Figure  6.  (This  figure 
as  well  as  Figures  7  and  8  are  not  drawn  to  scale;  they  merely  show  the  general  appearance  of  the  curves.) 

i  /  _2  2 

If,  however,  v  is  smaller  than  yl8  e  -  4c,  then  v  -  4p'  is  zero  in  two  points  which  lie  between 
y*  and  yS  and  to  different  sides  of  y  *  1/2.  Eetween  these  two  points  D  is  negative.  If  this  interval  excludes 
the  point  yj,  i.e.,  If  D(y^)  is  positive,  then  the  shape  of  the  curve  q'  =>  0  is  as  shown  in  Figure  7,  The  shape 
of  this  curve  is  further  modified  in  an  obvious  way  if  D(y^)  is  negative,  or  if  it  is  positive  but  y^  is  larger 
than  1/2.  In  Figure  7,  the  sign  of  q*  (y,  z)  is  indicated  by  +  and  -  signs.  In  the  +  regions  the  integral  curves 
cf  Equation  (7)  are  curved  upwards,  in  the  -  regions  they  are  curved  downwards.  Compare  this  diagram 
with  Figure  4  where  some  of  the  Integral  curves  are  shown. 

Figure  8  shows  the  types  of  curves  q‘  »  0  which  can  occur  in  the  case  c  »  0  for  various  values  of 
v.  The  dotted  line  corresponds  tov  *  4/e,  for  which  D  vanishes  only  at  y  .»  1/2  and  y  «  0. 

Let  us  now  calculate  the  derivative  of  q'  (y,  z)  with  respect  to  y  In  the  direction  of  the  integral 
curves  of  (7).  This  is 


q"  (y.  z)  -  -§y  q'  (y.  z)  +  ^  q’  (y>  z)  ■  q  (y. z)  ■ 

An  easy  transformation  shows  that 

q"(y.  *)•  )  q'  (y,  z)  . 

Now  p"(y)  is  positive  between  y  »  0  and  y  »  1/2,  and  negative  for  all  larger  values  of  y.  If  we  consider  a 
point  on  the  curve  q*  =  0,  then  q"  »  -p"/z  is  positive  for  y  <  1/2,  z  < 0,  and  for  y  >  1/2,  z  >  0;  q"  Is  negative 
for  y  <1/2,  z>0  and  for  y  >1/2,  z  < 0.  From  this  fact,  a  number  of  conclusions  can  be  drawn  concerning 
the  shape  of  the  integral  curves  of  (7).  For  instance,  consider  the  case  of  Figure  8  for  small  v.  The  Inte¬ 
gral  curves  of  Equation  (7)  have  positive  curvature  (q‘  >  0)  Inside  the  loop  which  is  formed  by  the  left-hand 
part  of  the  curve  q*  «0,  as  well  as  to  the  right  of  the  right-hand  part  of  that  curve.  If  an  integral  curve 
passes  through  a  point  insitls  the  loop,  then  the  entire  portion  of  the  integral  curve  to  the  left  of  that  point 
is  also  within  the  loop.  In  other  words,  an  Integral  curve  can  leave  the  loop  only  in  the  direction  of  increas¬ 
ing  y,  but  not  in  the  direction  of  decreasing  y.  For  in  a  point  where  the  Integral  curve  crosses  the  boundary 
of  the  loop,  q'  is  zero  and  q"  is  negative,  therefore  to  the  right  of  this  point  on  the  integral  curve,  q'  Is 
negative,  and  thus  the  points  on  the  right  are  outside  the  loop;  that  Is  to  say,  in  crossing  the  boundary  of 
the  loop  in  the  direction  of  increasing  y  along  an  integral  curve  of  Equation  (7),  we  are  leaving  the  loop, 
not  entering  it.  Similarly,  an  integral  curve  of  (7)  can  only  enter  but  not  leave,  in  the  direction  of  lncreas-  ' 
ing  y,  the  area  to  the  right  of  the  right-hand  part  of  the  curve  q'  -  0,  Analogous  conclusions  can  be  drawn 
about  the  various  portions  of  the  area  of  positive  q'  in  the  cases  of  Figure  0  and  Figure  7. 
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A  further  conclusion,  in  the  case  of  Integral  curves  going  through  points  Inside  the  loop  in  Figure  A, 
is  that  ail  such  curves  will  go  through  the  origin,  and  that  in  the  neighborhood  of  the  origin  their  curvature 
remains  positive.  Their  slope  at  the  origin  is  0.  There  are  infinitely  many  curves  of  this  kind,  since  they 
fill  the  entire  area  of  the  loop.  They  are  the  curves  of  "Type  1-A"  mentioned  In  "Solution  of  the  Transformed 
Equation  without  Heat  Loss,"  of  this  report.  A  similar  reasoning  applied  to  the  right-hand  portion  of  the 
area  where  q‘  is  positive,  leads  to  the  curves  of  "Type  III-C." 


SUMMARY 

Steady-state  burning,  as  defined  by  Equation  {5}  or  (6)  can  tatce  on  one  of  the  following  forms: 

a.  In  the  case  without  heat  loss  to  the  surrounding  space,  there  are  infinitely  many  solutions  In 
which  the  temperature  approaches  absolute  zero  asymptotically  as  x  tends  to  Infinity,  has  a  mrulrmirn  at 
some  finite  x  and  decreases  toward  zero  to  the  left;  one  of  these  Is  distinguished  by  the  fact  that  Its  maxi¬ 
mum  temperature  is  higher  than  that  of  any  other  solution  in  this  group.  There  are  also  infinitely  many 
solutions  (these,  however,  do  not  fulfill  the  boundary  conditions  of  steady -state  burning'  in  which  the  tem¬ 
perature  reaches  zero  for  a  finite  x  and  increases  beyond  all  bounds  as  x  decreases.  Between  these  two 
groups  of  solutions  there  Is  a  third  whose  character  depends  on  the  value  of  the  burning  velocity,  v.  If 
this  is  smaller  than  ,90280,  the  solutions  in  this  group  have  a  maximum  temperature  for  some  value  of  x, 
and  the  temperature  declines  toward  absolute  zero  to  both  sides  of  the  maximum.  If  v  Is  greater  than, 

or  equal  to,  .90280  then  the  temperature  In  this  Intermediate  group  of  solutions  approaches  Infinity  as  x 
decreases,  and  approaches  zero  as  x  grows  to  Infinity.  The  limiting  case  of  the  first  group  appears  to  be 
the  most  interesting  one,  and  several  examples  (for  various  values  of  v)  are  plotted  on  Chart  i3.  However, 
all  the  solutions  of  the  first  group  are  consistent  with  the  conditions  of  steady-state  burning, 

b.  In  the  case  with  heat  loss,  there  is  a  unique  solution  which  is  consistent  with  the  conditions  of 
steady-state  burning.  On  it,  the  temperature  reaches  a  maximum  for  some  finite  x,  declines  to  both  sides 
and  approaches  the  "initial  temperature"  yQ  (close  to  room  temperature)  asymptotically  as  x  grows  to  in¬ 
finity.  (For  declining  x,  the  temDerature  declines  to  absolute  zero;  but  this  part  of  the  curve  no  longer 
represents  burning  in  the  same  substance.)  Examples  of  such  solutions,  for  varying  values  of  flame  speed 
(v)  and  room  temperature  (y  )  are  shown  in  Chart  t3.  The  solutions  are  practically  Independent  of  the 
value  of  c,  the  constant  which  determines  the  amount  of  heat  loss,  as  long  as  c  is  reasonably  small  but 
different  from  zero. 

There  are  numerous  other  solutions  of  the  differential  equation  of  steady -state  burning  with  heat 
loss,  but  none  of  them  satisfies  the  initial  conditions  usually  associated  v:i*h  burning;  in  particular,  there 
Is  no  other  solution  which  starts  at  room  temperature.  There  are  Infinitely  many  solutions  which  approach 
the  ignition  temperature  y^  from  above  or  below  as  x  increases  to  infinity  {cirvss  of  Type  I-A  and  I-C, 
respectively);  others  are  similar  to  the  solutions  of  Type  II-C  and  ni-C  of  the  case  without  heat  loss  in 
that  the  temperature  increases  to  infinity  with  decreasing  x;  a  few  solutions  (shown  in  Diagram  4  near  the 
origin)  remain  at  low  temperatures  throughov.. 
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Chart  4  shows  all  types  of  solutions  (of  the  "transformed  equation")  for  the  case  In  which  the  para¬ 
meters  v,  c,  y  have  realistic  values;  In  particular,  that  v  Is  smaller  than  V,  that  y  Is  sufficiently  small 
a  s 

so  that  p(y)  has  three  real  roots,  and  that  c  is  small.  Other  types  of  solutions  occur  for  values  of  the  para- 

i  p 

tneter3  outside  these  limits.  For  Instance,  if  c  is  large  enough  relative  to  v,  so  that  p*(y^)>  ,  then  the 

transformed  equation  has  solutions  which  spiral  around  the  point  y^,  and  the  corresponding  solutions  of 
the  equation  of  burning  oscillate  about  the  ignition  temperature. 


t, 


Franz  L.  Alt 
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